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Localisation of Géométrie Anisotropy 
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Abstract 

The class of 2-D nonseparable geometrically anisotropic localisation operators is defined, containing highly 
~ anisotropic nearly unidirectional localisation operators, as well as isotropic localisation operators. A continuum of 

anisotropic operators between the extrêmes of near unidirectionality and isotropy are treated in a single class. The 
eigensystem of any given operator in this family is determined, thus specifying geometrically anisotropic optimally 

■ concentrated functions, and their degree of localisation. 

x> : 

Index Terms 

Localisation operator, wavelets, anisotropic and directional variation. 

f-H | I. Introduction 

^ r ■ i His correspondence introduces a new class of spatially anisotropic and nonseparable 2-D localisation operators, 
' _M_ namely the class of geometrically anisotropic localisation operators. The opération of 'localisation' in this 

■ context refers to limiting variation in a 2-D square integrable function to variation associated with a given set of 
^ | spatial locations and spatial frequencies, see also Daubechies [l][ch. 2]. By formally using a localisation operator, 

l— ~ 1 ' the energy concentration of a given function to an anisotropic and nonseparable 4-D région of 2-D space and 2-D 
— i spatial frequency, may be exactly quantified. Optimally localized functions can then be derived. 

The important set of tools which motivâtes the need for deriving optimally localised functions is that based on 
the local représentation of functions. A Fourier Transform (FT) represents a signal globally, Le. it décomposes a 
çnJ ■ signal in terms of modes présent over the full length of a signal. Given many signais exhibit transient features, it 
CN ", is necessary to represent the local properties of the signal. Examples of local décompositions include the Wavelet 
^ Transform (WT) and the windowed FT [2]. 

O ■ To be able to form a local décomposition of a given function, a set of well-localised décomposition functions 
! must be used, or the utility of the décomposition vanishes. The locality of a given décomposition function must 
be adjusted to the class of signal that will be analysed. By the construction procédure used for the WT [2], the 
locality of any member of the family of functions used for analysis, is in some sensé 'équivalent' to the locality 
of the mother wavelet. Thus it is sufficient to détermine appropriate mother wavelets for a given problem. The 
eigensystem of a localisation operator corresponds to a whole set of mutually orthogonal and optimally localized 
^ ■ functions, where the eigenvalues of the functions, measure their degree of localisation to a given localisation région 
[3], [1]. Thus eigenfunctions of localisation operators are suitable mother wavelets. 

By determining the eigensystem of the operators proposed in this correspondence, we obtain classes of mother 
wavelets suitable for analysis of non-stationary geometrically anisotropic fields. Thèse facts justify the construction 
and study of the proposed operators. 1-D localisation operators have already been the focus of considérable study, 
see for example the références in [1], [4], [5]. Consécutive truncations form a possible method of constructing 
localisation operators: such procédures treat the space and spatial frequency variables inhomogeneously. This results 
in an unequal compromise between concentration in space and in spatial frequency, a clearly undesirable feature. 

Also it is désirable to define operators localising functions to régions that do not contain zéro frequency, if we wish 
to build for example families of wavelets. In this case the relocation in frequency is not implemented by frequency 
shifts. Daubechies and Paul [4] defined a set of homogeneous localisation operators, that treated time and frequency 
variation on an equal footing. Thèse operators were used to dérive optimal 1-D mother wavelet functions, and the 1-D 
Morse wavelets were thus obtained. The Morse wavelets have been used to analyse geophysical, astrophysical, and 
médical time séries [6], [7], and have been extended to radial wavelets in 2-D [3]. The construction of localisation 
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operators in both cases started from forming a local décomposition of the observed function, discussed in section 
M 

We seek to develop operators for the study of géométrie anisotropy, formally defined in section HII-Ai Our 
motivation for studying this class is that it contains isotropic functions, and anisotropic functions that are nearly 
unidirectional, as well as the full continuum of structures in between isotropy and unidirectionality. Géométrie 
anisotropy is used for modelling covariance structures in geostatistics (see for example Christakos [8][p. 61]). Thus 
the optimally localised geometrically anisotropic mother wavelets will have a natural application area, in the study 
of non-stationary geometrically anisotropic random fields. 

In gênerai it is very hard to define operators so that their eigensystem may be determined [l][p. 41], and without 
a known eigensystem the operators lose most of their utility. We shall therefore define the anisotropic operator with 
great care, in a séries of steps, in section [TTT1 We start by defining the appropriate 'fiducial vector' from which 
to build the operator. Géométrie anisotropy is characterised by a transformation matrix, and thus a transformed 
distance metric is defined. The fiducial vector is picked to be a function marginaily well localized in space in terms 
of the transformed distance metric, and also well-concentrated in transformed frequency. From the fiducial vector a 
family of cohérent states must subsequently be defined. We construct the family of cohérent states in section IIII-BI 
using the transformation matrix of the géométrie anisotropy. In gênerai the family of cohérent states will not be 
constructed by the usual 2-D opérations of translating, scaling and rotating the fiducial vector [3], but rather by a 
set of transformed opérations. 

We subsequently need to détermine the localisation of an arbitrary élément of the family of cohérent states, as this 
in conjunction with the choice of région of intégration, will détermine the localisation of the operator, see section 
IIII-CI This enables the définition of operators with given localisation régions, see section IIII-DI We subsequently 
détermine the eigensystem of the proposed localisation operators, see section [ÏVj We dérive explicit forms for both 
eigenfunctions and eigenvalues. We give spécifie examples, for given transformation matrices, of functions in this 
class. 

Thèse developments combine to define new classes of 'optimally' localised nonseparable anisotropic functions. 
The developments span previously derived optimally localised functions, and form a cohesive framework for treating 
geometrically anisotropic 2-D localisation. 

II. Local Décompositions 

We shall construct the localisation operator starting from a local décomposition, similar to the 2-D Continuous WT 
(CWT), see for example Antoine et al. [2]. We define a family of 'cohérent states' f^(x) from the 'fiducial vector' 
[3], a single 2-D function w(x). w(x) is assumed to be marginally well-localised in space and spatial frequency. 
The index £ = [a, 9, h] T eief 4 régulâtes the localisation of the éléments of the family of cohérent states. For 
example, the usual WT corresponds to taking: u^(x) = VqTIqT^ {v} (x), where for a > 0, V a v{x) = a~ x v (x/a) 
is the dilation operator, for 9 G [0, 2tt), 1Zgv(-x) = v (r_0x) is the rotation operator, with rg as the rotation matrix 
[2], and for b G R 2 , T h v(x) = v(x — b) is the translation operator. 

Let u(x) have a FT given by V(u), in angular frequency uj. We do not refer to u(x) as a 'mother wavelet', 
as we intend to introduce a différent set of opérations to construct f^(x). We define the 'local coefficients' of 
v,g G L 2 {R 2 ) by: 

Wv&9) = j ^|(xMx) d 2 X = <^(x) )9 (x)> = j J V£(u>)G(u>) d 2 u. (1) 

Eqn. (OQ) represents a projection of g(x) into 'local contributions' w v (£;<?). The function <?(x) can be reconstructed 
from the local coefficients, if £ and u(x) satisfy a suitable set of constraints. For example if v(x) is an admissible 
mother wavelet, with admissibility constant C v and £ corresponds to the indexing denoting the translation, scaling 
and rotation opération, then we may reconstruct g(x) by: 

p(x) = — f vç(x)w v (Ç;g) dAç, dAç = a~ 3 da d 2 h dO. (2) 

The interpretability of eqns. £j]) and (O dépends on the locality of ug(x), in turn determined from the locality of 
v(x) and the choice of opérations denoted by Depending on the type of function g(x) that we are decomposing, 
différent families of cohérent states are suitable to use for the analysis of this function. To focus our interest on a 
spécial class of anisotropy, we shall now introduce the class of géométrie anisotropy. 
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Fig. 1. The n = eigenfunctions in space when a) P H = ([1 0] ; [0 0.15]) , b) P H = ([1 0] ; [0 0.5]) and c) P H = ([1 0] ; [0 1]), 
with 7) = (8, 3). The same functions are also plotted in frequency, in d), e) and f). 

III. CONSTRUCTING THE GEOMETRICALLY ANISOTROPIC OPERATOR 

A. Geometrically Anisotropic Functions 

Définition 3.1: A Geometrically Anisotropic Function 
A function g a (x) is said to exhibit géométrie anisotropy if for a fixed non-negative 2x2 symmetric matrix H and 
g r (-) a 1-D function, it takes the form: <7a(x) = g r (x T Hx) . 

An application of such functions is to model auto-covariance of random fîelds having géométrie anisotropy, see 
Christakos [8][p. 61]. For such random fields the auto-covariance of the field h(x.) at spatial locations x£R 2 and 
at y G R 2 takes the form bh (x, y) = cov {/i(x), h(y)} = b r ^(x — y) T H (x — y)J . An example of such functions 
is an isotropic covariance, given by H = I, where &h(xi,X2) is specified as a function of the Cartesian distance 
between xi and X2, but also for certain observed phenomena it is natural to model the auto-covariance between 
the field at two points as depending on a local affine transformation of the two variables [8], [9]. Note that the 
anisotropic extrême of this form corresponds to having H = Ii(e) = ([1 0] ; [0 £ 2 ])- The function becomes 
approximately unidirectional, as e — ► 0. 

We form the décomposition of H as H = P^Ph, and refer to Ph as the 'transformation matrix' of H. 
Unfortunately this spécification is not sufficiently constrained to uniquely détermine Ph, however as we in this 
correspondence only seek to détermine the optimal décomposition for a set H, any of the Ph will work. For H = I, 
we have P H = I whilst for H = Ii(e), P H (e) = ([1 0] ; [0 e]) . Let Q H = Ph - Define the 'transformed 
spatial variables' y and the 'transformed frequency variables' v by: y = Phx and v = Qh^- Then x = QhY 
and u) = Pgf. For convenience we also define y = ||y||, and v = \\u\\, to complément x = ||x|| (the modulus of 
the spatial coordinate) and oo = ||u>||. 

B. Defining the Geometrically Anisotropic Cohérent States 

We shall now construct a localisation operator that measures the localisation to régions that spatially decay in 
y 2 (x), and decay in frequency u(u), away from some set of frequencies. We first define the fiducial vector, very 
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much in analogue with choosing a mother wavelet, and from the fiducial vector define a family of cohérent states. 
Define as in Metikas and Olhede [3][eqn. 4], V^- D ^'^{u)) = ^—^Ç^-ujf 3 'e~^ for u > 0, and for u < 0. We 
use the radial cohérent states used in [3][eqn. 8]: 



00 uj- l / 2+ "i/ 2 



— ooJ — oo 



F (l-I3,A7)( a; ) e ^x d 2 aJ (3) 



Vf 



(/3.7), 



x) is marginally well localised, see [3]. Define the geometrically anisotropic fiducial vector by: 



vh(x) = y/\P^\v r (y), Vn(u) = |P H r 1/2 K(Qh^) = |P H 



,-1/2 



V r (u). 



(4) 



Let uq = arg w>0 maxw 



-1+7 IV 



(1-D,/J, 7 ), 



(lo)\ 2 . ujr' 7 ^(x) is decaying radially in space from x = 0, and is supported 



1D 



at angular frequencies such that u = u:q, where oj « ujq. Thus «h(x) is local in terms of x Hx and is associated 



with frequencies u such that | Q|j^ 

We now need to define the cohérent states from the fiducial vector defined in eqn. (01). We define the transformed 



.Thus if Ph is itself a 



rotation operator K 6 ,h by: 7^ iH M (x) = # (r_0 jH x) , r_ 0iH 

rotation, then Tv^ h corresponds to a normal rotation, otherwise TZq,h defines the act of rotation in the transformed 
space, i.e. in the y coordinates. 

We define the generalized spatial shift, also in the transformed space, for any function g(x) with FT G{u) by: 

1 



T^GH = G(«)e 



-jb T w 7 



(5) 



(2vr) 2 7_ c 

Thus if 7 = 1 this simplifies to the usual spatial shift, if by Qiib rather than by b. The generalized spatial shift is 
thus adapted to the transformed geometry. If 7 ^ 1 then the local effect on <?(x), local to any given wavenumber 
u = uq, is a shift in space by |Ph^o| 7 1 Qiïk. The magnitude of the shift dépends on |Ph^o| 7 1 • The act of 
a implementing a generalized spatial shift corresponds to shifting the fiducial vector in space by différent amounts 
at différent frequencies, in précise analogue with the radial [3] and 1-D [4] generalized spatial shifts. Introducing 
the generalized spatial shift enables us to treat a larger class of localisation régions. A family of cohérent states is 
then defined by: 

W x ) = ZV/ï^h^Vh (x) , W) = a 1/7 ^n(a 1/7 ([f-e.nr 1 ^)' 

1/7 



-j'u> T Q H b||Q H !*>| 



e -i^-V/7|p H |-l/2 



V r (a L 



and V&(u) = e-^+^^a^n^^nl 



7.V7 



r-0"), 

2/3 + -/-1 



(6) 



/(V|P H |r(r)). 



C. The Localisation of the Cohérent States 

Before expressing any function in terms of the family of cohérent states, and proving that a resolution of identity 
may be achieved, we shall détermine the localisation of the cohérent states. We need to détermine the spatial and 
spatial frequency locality of v^J (x) as a function of £, to be able to détermine the localisation of the operator. We 
firstly détermine the energy of the cohérent states (see eqn. (|22l) in the appendix), and noting that d 2 u; = |Ph| d 2 u, 
this then yields: 



(1>h(0 



(n 



09.7) ,,(0,7) 
H;£ ' H;£ 



2 1 r°° 

d 2 u = — uP 
^ Jo 



du. 



We then calculate the spatial and spatial frequency average of the function, just like Daubechies and Paul [4] [eqn. 2.10] 
. It is easiest to calculate the transformed coordinate averages, rather than the averages of x and u. We détermine 
that (see eqn. d23l ) in the appendix): 



<">h(0 



1 



<1)h(0 (2vr) 2 



v |Vh ; ^(^)| 2 d 2 u 



|Ph| 



7 2 /7 



(l)H(O 7T» 



V H (a 



1/7, 



3 u) 



d 2 v 
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Similarly we can détermine that the average spatial position of any member of the family of cohérent states is 
given (see eqn. (l24l ) in the appendix) by: 



/ \ tes 1 f°° f°° I t a i 2 1 ^- 2 ( 7 + l)r(r-l/7 + 2) 

{yi)n(î) = 777-777 / m Vh ; é(x) dx= 777-— 77 « /7 ^- 

<1>h(«) 7-00 r(r + i) 



We define C/f 3 ' 7 - 1 = 2 1 / T_2 (7 + l)r(r - I/7 + 2)/r(r + 1), to be the term multiplying a" 1+1 / 7 fy. The average 
spatial position (x)h(£) bas components (xj)h(£) = [Qh ({ui)h($) (2/2)h(0)]j • With this set of relationships 
we may note the locality of wj^' 7 ' (x) and use this to construct a suitable localisation région for a function exhibiting 
géométrie anisotropy. 



D. Defining the Localisation Operator & Resolution of Identity 

1) General Operator: For any g(x) G L (R) we define for any région icK + x [0, 2vr) x R 2 the localisation 
operator. Define firstly the région of space and spatial frequency given by: 

P H = {(x,u,): x=(x) h (É), |QÏH = Mh(É), (8) 
Then subsequently define the projection operator by: 

Ja 



G 



H 



(2vr) 2 
C H 



OO /"OO 



'.4 J — 00 J —00 
r>00 /"OO 



2tt 




.4 J —00 J —00 



a 2/7T/(A7) (a l/7 



M + x [0, 27r) x M 2 = _4. a n and note from Metikas and Olhede [3] [section HI.D] that assuming 
-27 d 2 v < 00, we have that p^G(wi) = G (P^i) = G (wi) . We thus achieve 

a 'resolution of identity'. This demonstrates that eqn. (© is a suitable localisation opération. If we include ail space 
when integrating then we retrieve the full function <?(x). Of course, to appropriately define a localised function, 
the operator needs to be specified very carefully. 

We start by defining the appropriate restriction of the local index £ for a fixed C > 1 by: 



We now let A 

= 2lF / iîfi 



A{C) = Ua,9,h): a 2 + 6 2 + 1 < 2aC, 8 G [0,2tt) , 6= ^ 2 + 6|| . 



(10) 



Note that A(C) does not dépend on H, the transformation of the distance metric, or the shape parameters (/?, 7). 
We shall now use the cohérent states defined by eqn. §6§ so that A(C) of équation ( fTÔl gets mapped onto a région 
of space and spatial frequency that is given by: 



yuE. 



(0.7) 



+ 1/7 - 



(0,7) 



< 



Sf 7)2 (C7 2 - 1) } , (11) 



with C^P n \ and C^' 7 '', constants whose values are given by eqns. (|7]) and ([8]). For simplicity we have also defined: 



E. 



(0,7) 



G 



(0,7) (7-1)' 



(0,7) 



E (P,7) _ ^0,7)7 _ 



(12) 



Thus the région ï>g' 7 ^(C) is specified by four différent constituent parts: the transformation Pjj, the shape 
parameters (/3, 7) and C, the hypervolume parameter. 



in D^' 7 ^(C) are therefore given by: 



<f ,7) / 



The local transformation Ph spécifies the local région '(C) in terms of the anisotropy between x\ and xi, as 



C - VC 2 - 1 

(0,7), 



cf 7 V 



We note from eqn. CEO) that C-VC 2 - 1 < a < C+ VC 2 - 1 and so th e maximum and minimum wave number 

> ^min 



c + vc 2 - 1 
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well as loi and L02, as the space is parameterised in y and u. The shape parameters (/?, 7) détermine the compromise 
between y and u, and C détermines the hypervolume of the région. We reparameterise the région T>^'' 7 \c) in y 



and v, as 'J (C). For a fixed value of C, the hypervolume of P^' 7J (C), is given by: 



(0.7), 



2? 



(/3,7) 
H 



(C) 



(i 2 x <i 2 tt; 



We change variables to a 

^(/3,7)Area(C) = 



a 



'(C) 



d 2 y d 2 u = 4tt 2 



y dy v du. 



1/1 > 



and 6 = y^&f+ftf 



(y3, 7X7-1) 



. Then: 



03.7) 
H 



(C) 



4vr 2 



c+ v / c^ï 



2aC -l-a= C (^) 6C (/3,7) C (A7) C (A7) da & 



.1-1/7 nl/7 



( 7 + l) 2 r 2 (r 



I + 2)r 2 (r + I + l) 



2 5 r 4 (r + 1)7 



2C\/C 2 - 1 + log 



7 

'C- v/C 2 " 



1 



c + vc 2 - 1 



We obtain the form starting from eqn. (1241 in the appendix. This dennes A(fl, 7) 
Area(C) 



( 7 +l) 2 r 2 (r-i+2)r 2 (r+i+l) 



and 



2 5 T 4 (r+l)7 

. Note that A(j3, 7) is only a function of (/?, 7) via r and so does 



not dépend on C, whilst Area(C) is only a function of C and does not change with ((3, 7). Also ^Area(C) > 
and so there is a 1-1 map between C and the area. 

Metikas and Olhede [3] [eqn. 22] defined a genuinely 2-D localisation operator, but were only able to détermine its 
approximate eigenfunctions. To be able to exactly dérive the eigenfunctions, they defined an operator only valid for 
radial functions, in [3] [eqn. 12]. Similarly, we intend to define an operator only valid for geometrically anisotropic 
functions. We assume that g(x) is a geometrically anisotropic function, and write it as <?a(x). Then in analogue 
with [3] we define the geometrically anisotropic cohérent state in terms of the simpler indexing of = (a, b) T by: 



■Cm 

,(l./3,7) 



(27T) 



V, 



(13) 
(14) 



As we may consider v^'ç'J' (x) as the cohérent state «ji'g ( x ) that ^as been averaged across directions in space 



over b constant (see [3]), we may consider v-j^'^J (x) local to the set of spatial positions 



x = Q H y: y=\(y)n(î)\, u> = Q H v: i/=(i/) H (0, ( 15 ) 

noting that |(y)H.(£)| 2 only dépends on a and b (see eqn. ©), whilst (z^)h(£) on ly dépends on a (see eqn. ^}). By 
only using (x) f° r e A-(C) with A-(C) = {(a, b) : a 2 + 6 2 + 1 < 2aC} we construct a function 

local in space and spatial frequency to 



XV >H (C) = {(x,u;): |Px| = |<y) H (£)| , |QÏH = Mh«), 



6 2 + 6 2 



(a, 6) G A(C)} 



In analogue with Metikas and Olhede [3] we then may define the localisation operator with (51, 52) H 
{2^)- 1 j G\{uj)G 2 {uj) v |P H | di/, by: 



P. 



(0.7) 
£>,,h(C) 



{G a } («) = Cn,A 



A r (C) ,6f * H a 3 



.(1./9.7) 



bdb. 



(16) 



Clearly by appropriate change of variables, we can note directly from Metikas and Olhede [3] [eqn. 14] that a 
resolution of identity may be achieved as long as Ch,a = {r — l)/2. 



STATISTICS SECTION TECHNICAL REPORT TR-07-01 



7 



IV. DETERMINING THE ElGENSYSTEM 

We now intend to demonstrate that the eigenfunctions of the proposed operators can in the instance of eqn. ([Toi 
be determined exactly, or in the instance of eqn. (© be determined approximately. 

Firstly we note that by eqn. (|9]) that determining the eigenfunctions of the operator defined from the Morse 
cohérent states for an arbitrary function corresponds to for a given Ph, (P,j) and C solving the set of eqns. 

p2?<?.7> (G) *(wi) = Atf(wi), or: 



I p h| 

S"i p poo poo 

= f f f «^^(^^r^Q^^e-^QS^IQS^r- 

(Ar Ja(o J-oo J-oo 



>A(C) J-oo J-oo 

^ 1 7)*( a i/7 r _ 0l/2 ) e ib T ^^->(p^ I/2 ) Sv 2 dA € . 
Of course we note that uj\ = P^i, and so we find that: 

^ 1 H ^ Ja(C)J -oo J-oo K { 

vV«>(a^r- g v 2 y' bTv * v 3~ 1 *(PÏ L v 2 ) d 2 u 2 dAç. (18) 

Let Z = /?— I and m = 7. Comparing eqn. (fT8T ) with [3][eqns. 7, 21], and denoting the eigenfunctions of the operator 
defïned by [3][eqns. 7, 21] as m (w), we obtain that the eigenfunctions of équation (fT8T ). denoted by V^^m h( u, )> 
are given by the same functional form, but where the argument has been adjusted to the géométrie anisotropy. 
^Ivim h( cj ) nave ^e same eigenvalues as the isotropic Morse wavelets, and we dénote them {A n;ri H(C)}. We may 
then note that the geometrically anisotropic eigensystem is given by: 

<1U,hH = *SU(QS«), A n;r)H (C) = X n , r (C). (19) 

Thus, whenever a set of eigenfunctions are determined for the operator of eqns. (7) and (21) in [3], thèse 
automatically correspond to eigenfunctions of the geometrically anisotropic localisation operator, once the argument 
has been adjusted. Determining the exact eigenfunctions of eqn. (fT8l ) is in gênerai not an analytically tractable 
problem. In analogue with Metikas and Olhede [3] we instead détermine the eigenfunctions of the geometrically 
anisotropic localisation operator: 

^{^Kx) = ÀVa(x) (20) 
V^ {c) m}M = C^f f t _^\ U ){^\g A )u^bdb 



Ar(C) 

C h ,a I ! ! « 



V-r(, /i; , 2 )-l/2+7/2 T/ (l-D,A7) ( „l/7, yi) COS W b ~ ^/ 4 ) 

/ O— 1-7TÏ / ' h 



A r (C) " \prhuJ{h 

*A (P^ 2 ) v(i-^) (fl i/^)«»^-f V2 dv2 d ^±. (21) 

Clearly comparing this with eqn. (13) of [3], we détermine that the eigensystem is given by: ^|,^ mH (w) = 

^11m(QH w ) an( ^ A n;rj H(C) = A n>r (C), where we note from Metikas and Olhede [3] that the eigenfunctions and 
eigenvalues take the form: 

= ^A n , l+xj%m v l e-» m Ll(2v™\ À n , r (C) = r^l^r-l) \T ^ " ^ ^ 

where c = (21 + 2)/m — 1 and L^(-) dénotes the generalized Laguerre polynomial [10][p. 783]. 

As an example of geometrical anisotropy for < e -C 1 we take: Ph(ê) = [1 0; e] , Qh(e) = [l 0; e -1 ] 
H(e) = [l 0; e 2 ] . We note that u = y/uf + e~ 2 u;|. We obtain increasing anisotropy as e approaches zéro. 
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In the isotropic limit we obtain the results of Metikas and Olhede [3] and with Ph = I, the isotropic Morse 
wavelets are retrieved [3]. Compare the range of possible localisation régions that may be found in this class: 
see Figure [TJ As we do not change the values of j3 and 7 thèse functions have the same localisation value for a 
fixed value of C. As we are for the highly anisotropic case extending the extent of the function in one direction, 
trivially, to conserve the area size, it is compressed in the corresponding canonical variable, as is apparent from the 
frequency domain plot. As the spatial extent extends in x%, the function squashes towards j% = 0. The proposed 
framework allows the quantification of the stretching explicitly, and the discussion of both e — » and e — ► 00. 



V. Conclusions 

This correspondence has defined géométrie anisotropic localisation, and derived the optimally local functions 
for geometricaily anisotropic and nonseparable régions of space and spatial frequency. A full set of eigenfunctions 
with associated eigenvalues have been determined for this class, useful for the study of non-stationary fields in 
geostatistics. 
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Appendix 

Renormalising the functions to unit energy using eqn. © we have: 
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We define cf ' 7) by Cf ' 7) = (ï/j)h(£)/(M 1/7_1 ) Furthermore the localisation région has size: 
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